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Motivation: Monte Carlo simulation and Sign problem 
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Exponential complexity to solve QFTs 


Quantum systems: the Hilbert space H and the Hamiltonian Å. 
Accept the lattice model of QFTs, then 


dimH = a. Volume | 


= Solving QFT completely is (expected to be) exponentially hard. 


Notice: We do not always need to know about all the states precisely. 
We are usually interested in 
Vacuum and several excitations, Thermal states, etc. 
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Path integral and Monte Carlo simulation 
To characterize the thermal state, we can use path integral, 
Z - try End = | DAexp(-5[4)). 
M 
The integration domain M (of an SU(N) gauge theory) is 
M= SU(N) Volume, 


We want a tool to evaluate Z without exponential complexity. 


Monte Carlo method: Consider the case e7314] > 0. 
Generate an ensemble {A;}; following e314, and evaluate 
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More about Monte Carlo simulation 


Do we really circumvent exponential complexity using MC method? 


- Typical values of |O(A;)| 
7 ; 
It indeed solves the exponential complexity for operators satisfying 


(O(A)) 
Typical values of |O(A;)| 


N 
1 
fad A, 
Error of N 3 O(A:) 


~ (B - Volume). 


It has been quite successful to understand Hadron structures, 
thermodynamics of finite-temperature QCD, etc. 


This is true only when e S4 > 0. 
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Sign problem and Exponential complexity 


To use Monte Carlo method when e14] 4 0, we generate the 
ensemble {A;}; following the phase-quenched distribution e12e(914). 


WA SID) a wD | O(A:)e GIA) 
Cai) MA Za + ye e-ilm(5[4]) 


(O(A)) = 


Since STA) œ 8- Volume, (eID), | = e72 Volume Af, Thus, 
Necessary N of configurations > e2? Volume As, 


Exponential complexity revives due to the sign problem. 


Question 
Can we make Af = 0 by inventing a clever technique? | 
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Example of sign problem: Finite-density QCD 


Partition function of QCD is 


Zaco(T.m) = | DA Det(P(A, pa) + m) exp (Sya). 


quark gluon 


Sign problem: Det(P(A, py) +m) #0 at uq £ 0. 
Neutron star 
ə Cold and dense nuclear matter 
@ 2Msun neutron star (2010) 


ə Gravitational-wave observations 


Binary NS merger (NASA) 
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Method: Path integral on Lefschetz thimbles 
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Lefschetz thimble for Airy integral 
Airy integral is given as 


3 
Ai(a) = [Zoi (= + ar) 


Complexify the integration variable: z = x + iy. 


4 


Integrand on R, and on A 
(a= 1) 
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Multiple integrals on Lefschetz thimbles 


Oscillatory integrals with many variables can be evaluated using the 


“steepest descent” cycles Js: (classical eom 5'(2,) = 0) 


f d"™x eS) — wam | dee"), 
M 


o o 


Unlike one-dimensional case, the steepest descent manifold is not 
uniquely defined. 
=> Use of the homology H„ (Mc, {eP < 1}) becomes quite 


essential: 
> ZIJ), 


XOZ [Ka]. 


K 


H,,(Mc, (ep Ret) < 1} 


Hp (Me (e0 YO < 1}) 


K 


[Pham, 1967; Kaminski, 1994; Howls, 1997] 
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Multiple integrals on Lefschetz thimbles 


Concrete construction Pick up a metric ds? = g,;dz' ® d2/, and 
consider the gradient flow: 


Jo are called Lefschetz thimbles, and Im[S] is constant on it: 


lim z(t) = co} l 
t--00 


Similarly, Ko = {z(0)|z(oo) = ze}. 
[Pham, 1967; Kaminski, 1994; Howls, 1997, Witten, arXiv:1001.2933, 1009.6032] 
[Christoforetti et al. (PRD(2012)), Fujii et al. (JHEP 1310), etc.] 


In = 4200) 
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Monte Carlo algorithm on Lefschetz thimbles 


Monte Carlo algorithms on thimble(s) has been developing: 


e Langevin on J, Cristoforetti, Di Renzo, Scorzato, 1205.3996 


ə Hybrid MC on J> Fujii, Honda, Kato, Kikukawa, Komatsu, 
Sano, 1309.4371 


ə Contraction algorithm Alexandru, Basar, Bedaque, 1510.03258 


ə Generalized thimble Alexandru, Basar, Bedaque, Ridgway, 
Warrington, 1512.08764 


These methods generate ensembles {z;}; on J, following e 


d”zi 
| NO) 
by d"zO(z) exp(-h715(2)) ~ N DDA 1 Jd” z;] i 
YA YA N LY 1 EULA 


N.B. satisfies resurgence for “nice” S (Berry, Howls, '91, Howls '97). 
To 


—S(z). 
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Applications 


Applications: Case studies of fermionic sign problem 
eo One-site Hubbard model (1509.07146, with Tomoya Hayata, Yoshimasa Hidaka) 


ə Multi-flavor massless QED% (1612.06529, with Motoi Tachibana) 
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Applications 


Case 1 One-site Hubbard model 
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AIET One-site Hubbard model 


Path integral for one-site Hubbard model 
We consider the (0 + 1)-dimensional fermion model, 


s= far (0/2 u= iole). 


The path-integral expression for the one-site Hubbard model 


(2 = 3 fy ard(r)): 


p f i 2 B2 
Za yf | do (1 + chip t+ ¥/2))? 9897/20, 
2mU R AB 


Fermion Det 


Integrand has complex phases causing the sign problem. 


p is an auxiliary field for the fermion number density: 


Yuya Tanizaki (RIKEN BNL) Fermionic sign problem June 28, 2017 @ ECT* 


16 / 27 


AJETET One-site Hubbard model 


Behaviors of number density and Lefschetz thimbles 
Number density n with exact result and one-thimble result: 


-0.5 0.5 1 15 
(YT, Hidaka, Hayata, 1509.07146)(cf. Monte Carlo with 1-thimble approx. gives a wrong result: 


> 
JU 
Fujii, Kamata, Kikukawa,1509.08176, 1509.09141, Alexandru, Basar, Bedaque,1510.03258.) 


Lefschetz thimbles with —0.5U < u < 1.5p: 


7 7 f T { x 
NA kh w aan N 
075 pabaya 
TT 
| \ 
Kilasi 
EY, + E + \ 
ose 27 cf a gi To t 
> GA ALE RON ATT eh ANAND PR PREP we df? a 
= CSET PR Ole tt ER AAA EDA 
e Se Che) Cee ee EER eee eee 7 
DO KA KAM TAL AT Ae TTP ERP PRE PE APA w/U + 1/2 
VEKAN THATS HD Ole TE PRR ee Pr PS 
ESPEN NETE EEA tafe rrr DA 
wo NN EA WA! Kalt t tot tiet tiea! tatt tr 
Yh ii TIN TANAN NAA ANT oe ATE 
\ | \ | | fi / 
075 UE 025 $ 025 05 Tis 
Re(z) 
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AIETE One-site Hubbard model 


Results 
Results for GU = 30: (YT, Hidaka, Hayata, 1509.07146) 
n 
ae Yi — 
sites 1-thimble | | 
--- 3-thimble 
— 5-thimble 


-0.5 0.5 “10 1.5 


u/U 
Necessary numbers of Lefschetz thimbles ~ BU /2r. 


(cf. Fujii, Kamata, Kikukawa,1509.08176, 1509.09141; Alexandru, Basar, Bedaque,1510.03258) 


Consequence 


In order to describe the step functions, we need interference of 
complex phases among different Lefschetz thimbles. 
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Cit a ales 
Semiclassical partition function 


Using complex classical solutions Zm, let us calculate the semiclassical 
partition function (YT, Hidaka, Hayata, 1509.07146): 


= 1 
Za = > e75m = e 50009, (« C KO 5) ad ji 


m=— CO 


This expression is a good approximation for —1/2 < u/U 8 3/2. 
ny Eny 


[U 
-1/2 0 1 Za 
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ATE One-site Hubbard model 


Sign problem after Lefschetz-thimble deformation 
This computation means that the sign problem exists after the 
Lefschetz-thimble deformation. 

Compute the reweighting factor: 


Fong, deze) Tes Br (p4N o?) 


|dz|e~Re(S@)) = = e7Re(Sm) g 05 (0, e-217/8U) 


At u = 0, for example, 


Is pa dze-S(:) U 


Consequence 


Lefschetz-thimble method reduces Af = 2U to Af = U/8 in the 
one-site Hubbard model. 
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Applications 


Case 2 Massless QED2 
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Meets GED 
Multi-flavor massless QED, 


2-dimensional U (1) gauge theories with Ny massless fermions: 


Ng 
Z= J DA e54 I DYDyexp | -X ii Pap, [Da - Ha] ba 
a=1 


Since the fermions are massless, the nonzero topological sectors do 
not appear: 

27 o, 27 1 

gt + z inde + *d@ + dr. 


o photon gauge 


A= 


o-dependence is computable using the anomaly equation, and does 
not have the sign problem. 


Yuya Tanizaki (RIKEN BNL) Fermionic sign problem June 28, 2017 @ ECT* 2/27 


ATE ntg Massless QED9 


Toron-field integral of multi-flavor massless QED» 
The toron-field integration becomes (7 = L/: temperature) 


27 
Z= i dhodh, exp |--Z Fto, m) A 
0 


where F is the fermion one-loop free energy (pw, = Lua/(27)), 


N 
1 ? T Lo 27 | 1 : 
F=N; (1-5) - EY fin (ear) 


621 n=l 


+ In (2 + gr AN) tn G 4 ee (nth —1+H4)4 a) 
Hin G + er rai) ae \. 


In the limit 7 — O, we can use the mean-field approximation with 
complex saddle points (1612.06529, with Motoi Tachibana) (cf.1504.02979, with 
Hiromichi Nishimura, Kouji Kashiwa). 
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Gradient flow for massless QED» 


All the relevant complex saddle points have real F: 


ZI NAY 
HIN Hs 


aaa 
bone 


Wat (Al 
BK NI 
NIN: 
zi YA Mi M N NI 

Figure: Gradient flow in the 'Re(hi)- rIm(h 0) plane at Ng = 3 
(1612.06529, with Motoi Tachibana) 


a a = = = nag 


Phase structure of multi-flavor massless QED» 


At the zero-temperature with finite L, 7 = 0, the first-order 
transition occurs in this model: 


Ng = Ng 
6h 
4 
2 


a 
1 3 5 
etre. 2 3 3 

(a) Fermion number at N; = 2 (b) Phase boundary at Ny = 3 
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Sign problem after Lefschetz-thimble deformation 


In the Ng-flavor massless QED,, there are only Ny relevant saddle 
points in the limit 7 - 0. 


Since the (complex) mean-field approx. is good in this limit after the 
Lefschetz-thimble deformation, one can show that 


Is. P Ah e-5(h) 


TG 7, [he RSW) ~1 


This implies that Af = 0. 


Consequence 
Lefschetz-thimble method solves the sign problem of massless QED». | 
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Summary 


ə The sign problem is reviewed from the viewpoint of exponential 
complexity. 


e Using Cauchy's theorem, one can deform the oscillatory integral 
into the sum of steepest descent integrals on Lefschetz thimbles 


e We consider two examples of the fermionic sign problem. 
In the one-site model, the exponential complexity is not solved, 
but Af is reduced. 
In massless QED», the exponential complexity is eliminated by 
using Lefschetz thimbles. 
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